1
Cơ sở lý thuyết và tính ổn định
MATH007Lesson 5
00:00
Trước khi tận dụng sức mạnh của các công cụ giải số như phương pháp Runge-Kutta bậc bốn hay công thức Adams-Moulton, chúng ta phải đặt ra một câu hỏi căn bản: Liệu nghiệm thực sự có tồn tại hay không, và nó có ổn định? Cơ sở lý thuyết của các bài toán giá trị ban đầu (IVPs) cung cấp dấu hiệu "xanh" về mặt toán học, đảm bảo rằng các phép xấp xỉ rời rạc hội tụ đến một thực tế vật lý có ý nghĩa thay vì nhiễu số học.

Nền tảng vững chắc: Tính liên tục Lipschitz

Để kiểm soát cách sai số lan truyền, chúng ta cần một hàm $f(t, y)$ không "nhảy" quá đột ngột. Điều này được hình thức hóa bởi Điều kiện Lipschitz.

Định nghĩa 5.1: Điều kiện Lipschitz

Một hàm $f(t, y)$ thỏa mãn điều kiện Lipschitz theo biến $y$ trên tập hợp $D \subset \mathbb{R}^2$ nếu tồn tại hằng số $L > 0$ sao cho:

$$|f(t, y_1) - f(t, y_2)| \le L|y_1 - y_2|$$

với mọi $(t, y_1), (t, y_2) \in D$. Hằng số $L$ này là "giới hạn tốc độ" cho sự thay đổi theo chiều đứng của hàm số.

Ví dụ 1: Phân tích các hằng số Lipschitz

Xét $f(t, y) = t|y|$ trên $D = \{(t, y) \mid 1 \le t \le 2, -3 \le y \le 4\}$. Theo định lý giá trị trung bình (hoặc tính chất của giá trị tuyệt đối):

$|f(t, y_1) - f(t, y_2)| = |t|y_1| - t|y_2|| = |t| \cdot ||y_1| - |y_2|| \le |t| \cdot |y_1 - y_2|$.

Vì giá trị lớn nhất của $t$ trong miền của chúng ta là 2, hằng số Lipschitz là $L=2$.

Tính toàn vẹn hình học miền

Chúng ta không thể giải một bài toán giá trị ban đầu trong một miền đầy lỗ hổng. Chúng ta cần Tính lồi.

Định nghĩa 5.2: Tập hợp lồi

Một tập hợp $D$ là lồi nếu với bất kỳ hai điểm $(t_1, y_1)$ và $(t_2, y_2)$ nào, đoạn thẳng xác định bởi:

$$((1 - \lambda)t_1 + \lambda t_2, (1 - \lambda)y_1 + \lambda y_2)$$

với $\lambda \in [0, 1]$ cũng nằm trong $D$. Điều này đảm bảo rằng không phần nào của đường đi nghiệm bị "rời khỏi" vùng tính toán hợp lệ.

Định lý tồn tại và duy nhất

Khi các điều kiện này phù hợp, chúng ta áp dụng Định lý 5.4: Nếu $f$ liên tục và thỏa mãn điều kiện Lipschitz trên một tập hợp lồi $D$, thì bài toán giá trị ban đầu $y' = f(t, y), y(a) = \alpha$ có một duy nhất nghiệm $y(t)$. Điều này chứng minh tính hợp lệ của các phương pháp đơn giản như Euler ($w_{i+1} = w_i + h f(t_i, w_i)$) hoặc phức tạp hơn như logic dự đoán-sửa lỗi:

$WC = w_{i-1} + \frac{h}{24}[9f(t_i, WP) + 19f(t_{i-1}, w_{i-1}) - 5f(t_{i-2}, w_{i-2}) + f(t_{i-3}, w_{i-3})]$.

🎯 Nguyên lý cốt lõi: Tính ổn định
Một bài toán là ổn định nếu một nghiệm duy nhất tồn tại và phụ thuộc liên tục vào dữ liệu ban đầu. Nếu hằng số Lipschitz $L$ cực kỳ lớn, bài toán trở nên "cứng". Trong các phương trình cứng, các phần chuyển tiếp suy giảm nhanh chóng, nhưng đạo hàm của chúng (độ lớn $c^n e^{-ct}$) thì không, điều này buộc phải sử dụng Thuật toán 5.8: Hình thang kết hợp với lặp Newton để duy trì tính ổn định.